Hall Resistivity and Dephasing in the Quantum Hall Insulator 
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The longstanding problem of the Hall resistivity p xy in the Hall insulator phase is addressed using 
four-lead Chalker-Coddington networks. Electron interaction effects are introduced via a finite 
dephasing length. In the quantum coherent regime, we find that p xy scales with the longitudinal 
resistivity p xx , and they both diverge exponentially with dephasing length. In the Ohmic limit, 
(dephasing length shorter than Hall puddles' size), p xy remains quantized and independent of p xx . 
This suggests a new experimental probe for dephasing processes. 
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The ideal Quantum Hall (QH) effect can be defined as 
the simultaneous quantization of Hall conductance and 
vanishing longitudinal conductance: 



G xy = ve 2 /h, G x 



0. 



(1) 



where v, the filling factor, is some integer or odd de- 
nominator fraction. Consider a perfect QH sample with 
four Ohmic leads. By Kirchoff's laws, it is easy to see 
that while the measured G xy ,G xx and R xx depend on 
the lead resistances, the Hall resistance is independent of 
them, and is precisely quantized at 



Rx V = h/ve 2 . 



(2) 



This robustness of R xy generalizes to larger circuits, as 
proven by Shimshoni and Auerbach [fjj for the two di- 
mensional Ohmic puddle network model. However, as 
we shall see below, the quantum interference between dif- 
ferent tunnel junctions (absent for the Ohmic transport) 
spoils the quantization of R xy [Q . 

Quantum interference also drives localization and re- 
lated T = QH to insulator transition ||||. This has 
been confirmed by explicit calculations of the localization 
length exponent, in quasiclassical approximation S and 
numerically ||,(7j on Chalker-Coddington (CC) networks. 

There is no consensus however, about interference ef- 
fects on the Hall resistivity; its value in the insulating 
phase is still controversial. Several theories expect it 
to be finite |||| P xy ° < 00 ; or even quantized [ ^0|JTT| 
Pxy° = h/e 2 ("semicircle law"). Experimental observa- 
tions also vary: some data can be fit by the Drude form 
fn^ , fl3f p xy oc B, while others see a weaker magnetic field 



(B) dependence, with p xy nearly quantized |1 - 

In contrast, Entin-Wohlman et. al. [[l5) pointed out 
that Px^ is not a self-averaging quantity in the insula- 
tor. Using a model of local hopping in external magnetic 
field, they concluded that p xy diverges in the limit T — > 0. 
Unfortunately, the applicability of this model to QH sys- 
tems is somewhat questionable. 
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FIG. 1. (a) Four-terminal L = 6 CC network. QH puddles 
are shaded. Edge currents with amplitudes Ap propagate 
along the arrows, (b) The incoming and outgoing amplitudes 
at the tunnel junction between puddles p and q are related by 
the scattering matrix 



The purpose of this paper is to determine the low tem- 
perature Hall resistivity in the magnetic field driven insu- 
lating phase. Differing from a previous calculation [pj, 
we directly compute the four-terminal resistances R a /3, 
a,[3 — x,y for random CC networks using the Landauer- 
Biittiker quantum transport theory. Distributions of R a p 
are obtained for different linear system sizes L and mag- 
netic field parameters 9, defined below. 

Our key result is that at large L, the Hall resistance 
R xy (8,L) of the CC network changes from the asymp- 
totically approaching quantized value h/ve 2 in the QH 
phase to an exponential divergence in the insulator [fl6|| . 
Second, effects of inelastic scattering due to electron- 
electron and electron-phonon interactions are introduced 
via a phenomenological dephasing length l^ (for 
convenience, we measure all physical length in the units 
of typical puddle size Zy). For l v shorter then the pud- 
dles' size (phase-incoherent limit) we map the chiral net- 
work to the Ohmic puddle network model , which has 
a quantized Hall resistance. In general, the macroscopic 
resistivity is approximately given by an Ohmic resistor 
network formed by elements of size L = l v with the com- 
puted distributions of R a p . The Ohmic network problem 
is solved using percolation theory following Ambegaokar, 
Halperin and Langer pi (AHL). This approach, justified 
by the wide distribution of R xx , yields the ensemble of 
elements which determine the Hall voltage. The macro- 
scopic Hall resistivity p xy is obtained as an average over 
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this ensemble. We show that the divergence of both Hall 
and longitudinal macroscopic resistivities in the insula- 
tor is consistent with one length scale £(0). We conclude 
with the suggestion that p xy {T) could serve as an exper- 
imental probe of the dephasing length. 

Four-terminal CC network. Consider a square CC net- 
work in Fig. [l], with corner puddles serving as external 
leads. Carriers, described at a given energy by complex 
amplitudes {A^}, where a,p are the edge and puddle 
labels respectively, can propagate only in one direction 
around the puddles, as indicated by arrows. The corre- 
sponding currents are just I p — \A p \ 2 . Edge amplitudes 
across each tunnel junction pq are related by a unitary 
scattering matrix s^ pq \ 

A^Ys^e^A"), s^=( COS °™ sin M ( 3 ) 



r'=p,q 



pq 



where 9 pq G [6 + 8,0 — 8], and edge phases <f>" £ [0, 2n) 
are independent random variables; we have chosen |7]] 
8 = to reduce the statistical errors. The magnetic 
field parameter 9 can be tuned across the QH-insulator 
transition at 9 C = 7r/4. In the QH regime 9 < 9 C the 
neighboring leads are connected by highly-transmitting 
edge states, while the insulating phase 9 > 9 C consists 
of weakly-connected puddles without global edge states. 
Near the transition, 9 — 9 C oc B — B c . Away from the 
transition, tunneling across saddle points yields asymp- 
totically tan(0) oc e -MB-B c ) Q_ 

The scattering relationships provide a set of linear 
equations for all A p which depends on the incoming am- 
plitudes A- n at the external leads i = 1-4. The global 
scattering matrix Sij[9 r , <j>"] relates the amplitudes at the 
external leads, A° ut — &ij so ^ e curren t 

transmission matrix — \ Sij \ 2 defines the probability 
of scattering from incoming channel i to outgoing chan- 
nel j. We assume that the external leads, also formed by 
an incompressible QH liquid, are in equilibrium, which 
implies that the chemical potentials of their edges are 

m 



tf = {h/ve)I?. 



(4) 



The resistance tensor is now given by a compactified ver- 
sion of the Biittiker-Landauer formula 

R a0 = (h/ve 2 ) A* a (l + T) P (P (1 - T) P)' 1 PA fj , (5) 
Al = (1,0,-1,0), A; = (0,1,0,-1), 

where the vectors A^ are transposed A Q , and the op- 
erator P projects out the zero eigenvector (1,1,1,1) of 
the matrix 1 — T (matrix S is unitary). It is easy to 
show that the field-antisymmetrized Hall resistance is 

Rxy — {Rxy Ryx) l^- 

Numerically, we calculated the LxL transfer matrix re- 
lating the amplitudes on L open wires at the left and right 



sides of the sample (see Fig. [I]). We used the transfer- 
matrix technique with intermediate orthogonalizations 
p2[ to reduce numerical errors for large systems. Then 
the closed edges were "linked" by eliminating the cor- 
responding pairs of incoming and outgoing amplitudes. 
This resulted in a reduced 4x4 transfer matrix, with 
which the full scattering matrix Sij was computed. The 
numerical accuracy was controlled by checking the uni- 
tarity of Sij , and by test runs at quadruple accuracy. 

To obtain the distributions of R xx and R xy for differ- 
ent values of 9 and the system size L, we repeated each 
calculation up to 10 6 times, taking (R xx ,R yy ) as two 
different members which have the same R xy . In the in- 
sulating phase 9 > 9 C , R xx has a very wide, resembling 
log-normal, distribution over several decades, as gener- 
ally expected in a localized insulator phase [^5|. The 
Hall resistance R xy is seen to have a much narrower dis- 
tribution, but with mean and variance also increasing ex- 
ponentially with the size of the system in the insulating 
phase. 

Phase-incoherent network model. In the opposite limit 
we consider a network where inelastic mechanisms com- 
pletely destroy quantum interference between different 
tunnel junctions. In this case different tunneling events 
happen independently, and the incoming and outgoing 
currents are related at each junction, 



s p l,\ 2 r3. 



(6) 



where the tunneling matrix elements s^ pq ^ are given by 
Eq. (||) . According to Eq. (||) , local chemical potentials at 
each edge are proportional to the corresponding currents. 
Thus, the dissipative and Hall voltages at each junction 
are related to the tunneling current I pq as 

V p f = R pq I pq , V p H q = sign(S) {h/ve 2 ) I pq , (7) 

where B pq = {h/ve 2 ) cot 2 9 pq is symmetric under rever- 
sal of magnetic field. These local relations reduce the 
network model to the Ohmic puddle network model of 
Ref. M, for which R xy {l v < ly) = {h/ve 2 ) for any real- 
ization of {9, (/)}. This model equally applies to integer 
and fractional QH regimes. Eqs. (Q) also imply that the 
phase-incoherent network model has an exact current- 
voltage duality M,O,0] , which interchanges the QH and 
insulating regions and simultaneously inverts the resis- 
tance R pq associated with each junction. This duality 
also inverts the macroscopic dissipative resistance of the 
system, determined, according to the AHL percolation 
argument pa], by the median of the distribution of R pq . 

To describe the regime of intermediate dephasing 
length, ly < l v < oo, where ly is a typical puddle size, 
we use finite-size four-lead CC networks. In the simu- 
lation, the phase breaking occurs only at the leads, and 
we define L — l v . Square blocks of linear sizes l v or 
larger connect as Ohmic resistors, with the resistances 
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Rap{L — lip) chosen from the numerically determined 
distribution. Ignoring relatively small Hall voltages at 
this stage, we notice that R xx is exponentially widely 
distributed. So, we can again follow AHL and insert 
all resistors in increasing order of R xx until the percola- 
tion threshold is reached. Here we assume a percolation 
threshold of 50%, applicable for self-dual lattices. Thus 
the last resistors to connect the percolating cluster (PC) 
have the median value Rf x d - Since all higher resistors 
R xx > R' xx d transport negligibly little current, they can 
be discarded. By AHL, and subsequent numerical confir- 
mations Q , the macroscopic resistivity of the network 
is simply p xx = i?™ ed (^). 

The Hall resistivity p xy can now be determined as the 
weighted average over the PC. Since numerically, \R xy \ <C 
R xx , and we draw no current from the Hall leads, it is 
safe to assume that the Hall voltage fluctuations across 
the PC are averaged out by secondary local currents, and 
the macroscopic Hall resistivity can be estimated by the 
average p x ™ ed = (R xy ) R xx <R™<> d ■ As an upper bound, 
we also calculate the average Hall resistance of the necks 
of the PC given by p™ ed = (R xy ) j^^^d , which gives 
similar results as shown below. 
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FIG. 2. QH to insulator transition in finite size scaling. 
Solid lines are the best fits of Eq. (^), and dashed lines are 
guide to the eye. Inset: The resulting inverse correlation 
lengths are consistent with the critical exponent v = 7/3 

Results. In the localized regime, 8 > 9 C) both p xx |?|] 
and p xy diverge exponentially at large l v . Both resistiv- 
ities fit very well [see Fig. |] for the fit of p xy (L)] to an 
exponential scaling function with finite-size correction 

p(L)~A e Liexp[L/Z(6)], (8) 

where 7 is a 0-independent exponent determined by the 
geometry of the system, and is the localization 

length (both quantities are defined separately for p xx and 
p xy ). As illustrated in the inset in Fig. |[ the divergence 
of both correlation lengths as 8 — > 8 C is consistent with 
£(#) ~ (8 C — 8)~ 7/>3 . This form for t; xx is in agreement 
with previous studies ||f7| . The fact that the correlation 
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FIG. 3. Asymptotic correlation between p xx and p xy deter- 
mined by two averaging procedures (see discussion in text). 
k is the slope on the log- log plot. 

length for p xy also diverges with the same exponent im- 
plies that the transition is actually governed by a single 
length scale £(#), and asymptotically p xy ~ (p xx ) k — ► 00 
with k ~ 0.32 — 0.35 (Fig. ||). This divergence, one of 
the primary results of this paper, contradicts the expec- 
tations p[-pl| of finite Hall resistivity in the insulator 
at l v — > 00. Our results agree with the conclusion of 
Ref. although here we have used a different model, 
applicable for transport in quantizing magnetic fields. 

Fig. U also shows p xy in the metallic QH phase for 
8 > 8 C . In this phase, as expected, p xy {L) approaches 
the quantized value at large system sizes. We did not at- 
tempt to determine the corresponding correlation length 
since our geometry has narrow leads (see Fig. [l]). We did, 
however, check that our results are not limited to systems 
with narrow leads by making a limited set of runs for a 
special self-dual network geometry, which remains iden- 
tical under the interchange of QH and insulating regions 
and the replacement 8 — > 28 c — 8. 

We also attempted to suppress the quantum interfer- 
ence by calculating the ensemble averaged (Ty), which 
is formally equivalent to considering the same non- 
interacting system at very high temperatures. This aver- 
aging resulted in both Hall and longitudinal resistance of 
insulating phase much smaller then the average quantum 
values. However, the specific values of these resistances 
differed significantly from the values obtained numeri- 
cally for the phase- incoherent networks of identical geom- 
etry; particularly, the Hall resistance was not quantized, 
with the offset increasing into the insulating phase. This 
demonstrates that at T > 0, without inelastic scattering 
quantum interference cannot be completely suppressed. 

In the fractional Hall effect regime, electron-electron 
interaction effects are primarily threefold: (i) Stabiliza- 
tion of fractional v < 1 QH phases in the puddles, 
(ii) Renormalization of inter-puddle electron tunneling 
rates, and (iii) Dephasing of charge carriers by inelas- 
tic scattering. At low temperatures the second effect is 
expected to be finite, since infrared divergence of the tun- 
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neling amplitudes is cut off by the finite size of the pud- 
dles. In this regime, the effective edge state transport 
theory for the fractional and integer cases is identical M , 
up to factors of v in the chemical potential relations (3). 
Thus it is legitimate to use the CC model as a model of 
coherent quantum transport on edges of fractional QH 
puddles. However the dephasing mechanism and deter- 
mination of lip remains an interesting open problem, for 
both integer and fractional cases. 

Quantum critical point and dephasing. In a vicinity of 
the true quantum critical point, the effective dephasing 
length should diverge at small temperatures as 

-l/z 



(9) 



Experimentally, resistivity and non linear resistivity data 
at transitions between plateaux have been collapsed onto 
universal curves using z « 1.0 (and independently deter- 
mined ~ 2.4) p3| , p^| , p4] | . Based on our results, at a 
quantum critical point, and in the absence of additional 
phase-breaking mechanisms, one would also expect the 
Hall resistivity to diverge as 



Pxy ~ (hjve 2 ) exp \kT~ 1/z {B - B, 



\7/3 



(10) 



Experiments, however, have reported a constant or 
weakly independent Hall resistivity on the insulating 
side of the transition. For similar samples, resistivity 
saturation at low temperatures has been reported [ p5[ . 
Evidently, both effects are inconsistent with a true zero- 
temperature quantum Hall to insulator transition, char- 
acterized by a diverging dephasing length (|{|). 

We conclude that a nearly quantized Hall resistivity 
indicates a strongly dephased regime where, i.e., l v < ly- 
In contrast, the longitudinal resistivity in this regime is 
expected to diverge at large fields as p9| , 

Pxx{B) ~ (h/e 2 ) exp \v {l v /lf (B/B c - 1)] , (11) 

where I is the Landau length. This expression allows an 
independent estimation of ly, and an upper bound on 
the dephasing length in the limit of small temperatures. 
Experimentally, the interplay between ly , the scale of the 
long-range potential fluctuations, and the phase breaking 
length l v was noticed in Ref. P6[ . 

It is not yet clear what mechanism can explain zero 
temperature finite resistivity in disordered QH systems, 
or why dephasing seems to be more pronounced in some 
particular experiments. One possible source might be a 
coupling of the edge excitations to nearby domains of 
compressible v = 1/2 phase. 
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